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I. INTRODUCTION 



When studying tunnelling in non-separable systems 
with more than one degree of freedom, one immediately 
encounters difficulties which generically can be traced 
back to the absence of enough constants of motion. Even 
in the very peculiar case of integrable systems, where 
continuous symmetries provide as many constants of mo- 
tion as degrees of freedom, as soon as separability is lost, 
the analysis of tunnelling is not a simple generalization 
of what occurs in one-dimensional (ID) autonomous sys- 
tems. The later case is detailed in textbooks on quantum 
physics (see for instance [[!]) and it has even been possi- 
ble to give a comprehensive analytical treatment in term 
of complex solutions of the Hamilton equations [|| . How- 
ever, it is not until the mid-eighties that a satisfactory 
quantitative approach has been proposed |^] for tun- 
nelling in non-separable integrable systems involving a 
larger number of dimensions. Moreover, integrability is 
a property of higher dimensional systems which is not 
generic. The coupling between several internal degrees 
of freedom as well as the coupling to an external source 
usually destroys some global constants of motion. With 
such a lack of constraints on the dynamics, the classical 
motion in phase space may become chaotic: it may ex- 
plore volumes with higher dimensionality and therefore 
becomes exponentially sensitive on initial conditions. It 
is not surprising that these deep qualitative differences 
between an integrable regime and a chaotic one appear 
at the quantum level too. Some of the properties of a 
quantum system do change when constants of motion 
are broken and, indeed, it is the very object of quan- 
tum chaos to study the signatures of classical chaos at 
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the quantum level (see for instance || to see how rich, 
vivid and successful this domain is). 

We define tunnelling as a quantum process which is 
classically forbidden to real solutions of classical equa- 
tions of motion. In this paper, we consider Hamiltonian 
systems only and study how the non-dissipative break- 
down of continuous symmetries affects tunnelling. We 
do not consider how tunnelling is modified by dissipation 
and decoherence of the quantum wave. Of course, this 
requires a great care in real experiments where making 
dissipation negligible is always a hard task. This is one 
of the main reasons why very few real experiments have 
been done on these questions, though they would defi- 
nitely help to understand tunnelling in the presence of 
chaos (as far as we know, the only experiments explicitly 
made on chaotic tunnelling in the XX th century are those 
presented in || with electromagnetic microwaves instead 
of quantum waves). 

During the last fifteen years however, theoretical and 
numerical investigations on autonomous 2D and time- 
dependent ID Hamiltonians systems highlighted some 
mechanisms ||, |i[ 0, ||, |[ [u| [ll], |l2) and much insight has 
been gained on the influence of non-separable dynamics. 
Experimental evidence of such mechanisms would be of 
great interest especially in the light of the subtle interplay 
between interferences and disorder. These phenomena lie 
in the general context of wave transport in complex media 
where the role of disorder is played by the (deterministic) 
chaotic dynamics instead of having a statistical random 
origin. Of course, other important motivations can be 
found in the numerous domains where tunnelling plays a 
crucial role as a fundamental quantum process like ion- 
ization absorption, nuclear radioactivity, molecular 
collisions, mesoscopic physics etc. More speculatively, 
studies on tunnelling in high dimensional Hamiltonian 
systems should provide us with a natural extension of 
the instanton techniques which deal with quantum field 
theories reducible to effective ID autonomous Lagrangian 
systems. 
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In 2001, it has been shown both theoretically IQ 
and experimentally jl5], [h| 0, |l8) that atom cooling- 
techniques (and possibly molecular physics as well, 
where formally similar systems have been extensively 
studied @ |§ ||) provide systems which fulfill all 
the severe requirements to study tunnelling in the pres- 
ence of classical Hamiltonian chaos: accurate manipula- 
tion of internal and external degrees of freedom, precise 
control of dissipation and decoherence and on the prepa- 
ration/detection set up. For a brief account intended for 
a large audience see [|4, |25|, ^6|). The present article 
has the ambition to mark out the future challenging ex- 
periments, that remain to be done for reasons which will 
hopefully appear clear in the following. 

This paper is organised as follows. In section || we 
give a general and informal overview. In section pTj| , 
we quickly recall the main theoretical apparatus that is 
needed in the following. We implicitly refer to |l4| for 
details and demonstrations. In section IV, we comment 
the results of and . In section [V , we show in this 
context, with the help of numerical experiments, the very 
precise form taken by the phenomenon known as chaos 
assisted tunnelling. We explain why it has not been ob- 
served yet with real atoms and propose how to actually 
bring it to the fore. Before the concluding remarks in 



section VII, we give in section ^ some more numeri- 
cal results which illustrate how subtle the signatures of 
chaotic tunnelling can be. 



II. CHAOTIC TUNNELLING 
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FIG. 1: A generalization of the paradigmatic double well po- 
tential (a) is to consider tunnelling between stable islands that 
are related by any discrete symmetry in phase space. Case 
(b) corresponds to Hamiltonian (^) with 8 — 1 and 7 = 0.018. 
Here, the time reversal symmetry plays the role of parity in 
case (a). 



The simplest situation to illustrate tunnelling is proba- 
bly the case of a particle placed in a ID time- independent 
symmetric double- well potential (see Fig. [j] (a)). Starting 
in one well with an energy that is below the maximum 
of the potential, a quantum particle can jump into the 
other well with a non zero probability though it is a for- 
bidden classical process. In addition to the classical time 
scale r given by the oscillating period inside one well, 
we therefore have a longer time scale, the tunnelling pe- 
riod T > t of the oscillations between the wells. In 
the eigenenergy spectrum, tunnelling appears as a quasi- 
degeneracy of the odd and even-symmetry states whose 
energies are both of the order of H/r but differ by an 
exponential energy splitting 

A, = ^ ~ e-*'* (1) 

where A is a ^.-independent typical action and can be 
interpreted in terms of a unique complex classical trajec- 
tory under the barrier ||, |27| . 

In the following we generalize this elementary situa- 
tion in two ways. First, unlike the parity in the pre- 
vious example, we can deal with a symmetry which is 
not necessarily either a spatial one or a two-fold one. In 
other words, we can have any discrete symmetry group 
acting on the whole phase space as well as any iV-fold 



symmetry which lead to bunches of iV-uplets in the en- 
ergy spectrum (or bands if N ^> 1). In the following 
we keep N = 2 since we have a two-fold symmetry T 
actually playing the role of parity (see Fig. |l| (b)) and 
being somehow decoupled from the other discrete sym- 
metries. The classical structure in phase space is globally 
invariant under T and the quantum eigenstates can be 
classified according to their symmetric or antisymmetric 
character under the unitary transformation which repre- 
sents T in the Hilbert space of states. Because T acts in 
phase space, it is usually more complicated than a pure 
spatial transformation. Thus, the two regions of phase 
space connected by quantum tunnelling, but classically 
not connected, are in general not separated by a simple 
potential barrier, but by a more complicated dynami- 
cal barrier. In such a case, tunnelling is emphasized to 
be called "dynamical tunnelling" as suggested by Davis 
and Heller in ||^]. It often happens that the classically 
unconnected region are associated to the same region of 
configuration space, with different momenta. A simple 
study of the density probability in configuration space is 
then unsufficient to characterize dynamical tunnelling; an 
analysis of the density probablility in momentum space 
is required. 

The second kind of generalization leads to much more 
puzzling questions. When dealing with systems with sev- 
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eral degrees of freedom or, equivalently, if an external 
time dependence exists, classical trajectories generically 
lose their regular behavior, cannot analytically be com- 
puted and are organized in a fractal hierarchy that is 
described by the KAM perturbative scenario. Recently, 
important progress has been achieved in the understand- 
ing of the continuation of these intricate structures in 
complex phase space and their role at the quantum level, 
see 1 29 , Q and especially Q . Anyway, we are therefore 
led to the following typical quantum chaos question: if 
one is able to create two symmetric stable islands sep- 
arated in classical phase space by a chaotic sea whose 
volume is under control (see Fig. |^), what is the effect 
of this sea on the (dynamical) tunnelling between the 
islands? 
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FIG. 2: Poincare surfaces of section correspond to Hamilto- 
nian (^) with 9 = 1.724137. Two stable islands in the vicinity 
of the origin are created by a pitchfork bifurcation at 7 ~ 0.56. 
Above this value, chaotic motion progressively invades phase 
space in between the two stable islands. At 7 > y c ~ 0.625, 
the latter are no longer connected by regular trajectories. The 
experimental configuration used in the nist experiments cor- 
responds to 7 = 0.96 just before the islands disappear in a 
bifurcation cascade at 7 ~ 0.97. 

The "dual" situation where chaos is created inside the 
wells while the dynamical barrier is kept regular has 
been introduced and studied theoretically and numeri- 
cally in jllj] . For a better understanding of what occurs 
in the energy spectrum when regular wells are separated 
by a chaotic sea, it has been proposed in Q to slightly 
break the tunnelling symmetry. Nevertheless, in the fol- 
lowing of the present paper, it must be kept in mind that 



a discrete symmetry will be always maintained exactly. 
At last, a third kind of generalization where the Hamil- 
tonian character is destroyed by introducing dissipation 
and/or coupling to a thermal bath, is beyond the scope 
of this work l|33fl . 



III. EFFECTIVE HAMILTONIAN SYSTEM 

Following @ [TJ, ^ (see also 1341 in a different con- 
text) we deal with an effective ID time-dependent system 
whose Hamiltonian is 



H(p, q; t) = — - 7 



+ cos t) cos q 



(2) 



in dimensionless units. 7 and 9 are two classical real pa- 
rameters that can be modified in real experiments. In 
addition, there is also one parameter, namely ft e ff, which 
fixes the quantum scale and is defined by the usual rela- 
tion between canonical operators: [q,p] — ih c g. It turns 



out that h e g is not constant any longer (see section IV be- 
low). It can be also experimentally varied via the rescal- 
ing factor that is needed in the canonical commutation 
relation in order to work in dimensionless units used to 
write (§). 

The time dependence breaks the conservation of en- 
ergy and therefore may generate chaos. In order to deal 
with such a Hamiltonian, it is crucial to remark that it 
has both a spatial and a temporal periodicity. The lat- 
ter implies that the Floquet theorem can be used, which 
states that the Hilbert space is spanned by an orthonor- 
mal eigenbasis of the evolution operator over one period. 
The corresponding eigenvalues of this unitary operator 
are distributed on the unit circle and therefore are la- 
belled by their phase, which is conveniently written as 
exp (i27re/7ieff), where 2-7T stands for the period of the 
modulation and e can be interpreted as a quasi-energy, 
a generalization of the notion of energy level for a time- 
periodic system. 

The spatial periodicity of the Hamiltonian is also ex- 
tremely important, as it makes it possible to split the 
Hilbert space into independant components, each com- 
ponent being characterized by the so-called Bloch vector 
k in the ] — 0.5, 0.5] range: under translation of 2ir along 
q, the quasi-energy eigenstates are just multiplied by the 
phase factor exp (i27rfc). One is thus reduced to solve the 
Floquct-Schrodinger equation in an elementary spatial 
cell with boundary conditions depending on k. One thus 
generates - for a fixed value of k - a discrete quasi-energy 
spectrum ej(fc). When the full range of k values is consid- 
ered, one obtains the familiar (quasi-) energy bands [[35). 

There is an additional discrete symmetry which can be 
used. The Hamiltonian (||) is invariant under the time- 
reversal symmetry (q,p,t) — > (q, —p, —t). In the classical 
surfaces of section, this implies a symmetry with respect 
to the q axis. In situations like the one in fig. [Tfs, this 
implies the existence of pairs of symmetric classically un- 
connected tori, i.e. a situation where tunnelling could be 
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observed. In the quantum world, the situation is slightly 
more complicated, because this symmetry connects the 
k subspace to the —k subspace. In the particular case 
k = (k = 0.5 could also be used), this implies that 
the Floquet eigenstates can be split in two subclasses of 
states which are either even or odd under the symmetry 
operation. The splitting between a doublet of even and 
odd states, Ae„ = |e+(0) — e~(0)| will be a measure of 
tunnelling. 

We will extensively use the Husimi representation of a 
quantum state pq |. Such a representation associates to 
each quantum state \ip) a phase space function ip H (p,q) 
(where p and q are real numbers) defined by 

r(p,q) = \(zm\ 2 (3) 

where \z) is the coherent state corresponding to the com- 
plex number z = (q + ip) / \f2h c g. Because \z) is a mini- 
mal Gaussian wave packet with average momentum p and 
average position q, the Husimi function tp u (p, q) contains 
some information about the degree of localization of \ip) 
in phase space, and makes it possible to associate quan- 
tum states with classical phase space structures. 



be the superposition of tunnelling oscillations (with dif- 
ferent frequencies) for various sets of parameters. This 
will at best - if the dispersion of the parameter values 
is reasonnably small - blur the oscillations at long time 
and at worst will completely destroy the signature of dy- 
namical tunnelling. It is experimentally rather simple to 
keep an accurate time-periodicity of the driving signal, 
i.e. to fix h c g. Similarly, the balance between the con- 
stant and the oscillatory term, hence the parameter 
is easily controlled. The 7 parameter is proportional to 
the laser intensity and may thus slightly vary across the 
atomic cloud (because of the transverse structure of the 
laser beams). The most difficult part is to control that 
a single Bloch angle k is excited. Indeed, this requires a 
phase coherence of the initial wavefunction over a large 
number of laser wavelengths, which is extremely difficult 
to achieve experimentally fl37fl , as will be shown in the 
following. In any case, the inhomogeneous broadening of 
the experimental signal because of the dispersion in k will 
be responsible for a decay of the tunnelling oscillations. 

A. NIST experiments @ 



IV. EXPERIMENTS WITH COLD ATOMS 

Under some severe conditions which constrain the ex- 
periments, Hamiltonian (^|) can be obtained as an ef- 
fective dimensionless Hamiltonian for cold neutral inde- 
pendent atoms of mass M interactin g w ith two coun- 
terpropagating laser beams [l^, [16|| . These two 
beams have two slightly different frequencies at lol+Slo/2 
and ujl — dw/2. The longitudinal coordinate x and 
the dimensionless q are related by q = 2k^x where 
fci = lul/c. The rescaling of the momentum is given 
by P — {2kL/MSu>)p x . 7 and 9 are fixed by the inten- 
sity of the lasers and the detuning of the laser frequen- 
cies with respect to the atomic resonance. The dimen- 
sionless time t is taken in Slu" 1 units and the expression 
of the effective Planck constant is h c ff = 8lor/5lo where 
lur = hk\j2M. Since in fine we want to measure ex- 
ponentially small tunnelling splittings Ae, it requires to 
maintain these conditions for a time at least larger than 
h c ff/Ae. Moreover, a very accurate control of the prepa- 
ration of the initial state and of the analysis of the final 
state is compulsory. 

As shown above, because of the temporal and spatial 
periodicity of the Hamiltonian, observing the standard 
signature of tunnelling, that is an oscillation of a quan- 
tum state between two classically unconnected regions 
of phase space requires that a single doublet of Floquet- 
Bloch eigenstates is initially populated, with well defined 
values of the parameters (7, 9, ft e ff) but also with a well 
defined value of the Bloch angle k. If more than a sin- 
gle doublet is populated, additional frequencies (related 
to energy differences between the various populated Flo- 
quet states) will appear in the temporal evolution. If 
any parameter is not fixed, the experimental signal will 



In the nist experiments, the two stable symmetric is- 
lands are chosen quite close in phase space in order to deal 
with not too small splittings. Another crucial point of 
this experiment is that the classical motion of the islands 
over one period, unlike those in ]l4| and 16 , always re- 
mains trapped in one spatial elementary cell of length 2n. 
The quantum states localized in these islands are conse- 
quently only weakly sensitive on the boundary conditions 
which are governed by the Bloch angle. In other words, 
the tunnelling period will be only weakly dependent on 
the Bloch angle k, which implies that the unavoidable 
broadening over k will not spoil too much the signa- 
ture of tunnelling. This is a major improvement over 
the tunnelling described in JTij and Jig ], where a very 
narrow band of Bloch angle is required to observe clear 
tunnelling oscillations. Moreover, the atoms involved in 
the tunnelling process stay longer in the region where the 
laser intensities are uniform. 

Indeed, as proposed in p8[ , [ p"7[ chap. 4 and 5] the 
two stable symmetric islands are created from a pitch- 
fork bifurcation of the fixed point at (p 7 q) = (0,0). To 
visualize it (see Fig. ^| (a) and (b)), we extract a one- 
parameter sequence by varying 7 while 9 is fixed to 
the experimentally chosen value in |l5|, i.e. 9 = 1.7. 
When 7 is increased, the pairs of symmetric tori appear 
at 7 ~ 0.56. At the center of each set ot tori, there is a 
periodic orbit. Over one period of the driving, the peri- 
odic orbit is essentially a rotation over the fixed point at 
(p, q) = (0, 0), which explains that the whole structure re- 
mains trapped in a single spatial cell. For 0.56 < 7 < j c , 
the tori remain nested in one connected stable island. 
At 7 = 7 C ~ 0.625, a chaotic sea separates the symmet- 
ric islands which shrink and move away from the central 
point before being dissolved through a cascade of bifur- 
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cations starting at 7 ~ 0.97. 

In one series of experiments, 7 ~ 0.96 and h c s ~ 0.8, 
the atoms are prepared in one island and their average 
momentum (p) is measured stroboscopically every mod- 
ulation period (— 2ir in dimensionless units). Since, in 
phase space, the islands rotate about the origin with the 
same period, if no tunnelling occurred no variation in (p) 
would be noticeable. In fact, starting the measurement 
sequence when (p) has its maximum value, oscillations 
are observed which illustrate the back and forth motion 
of the atoms between the islands due to dynamical tun- 
nelling. The tunnelling period T is about 10 modulation 
periods in this case (200 fjs). This is in perfect agreement 
with the quasi-energy splitting obtained numerically for 
the two Floquet eigenstates having the largest Husimi 
functions inside the islands. 

It is worth noting that the NIST group uses a Bose- 
Einstein condensate as a preliminary step for prepar- 
ing atoms in well defined quantum states, especially for 
achieving a large coherence length for the wavefunction, 
i.e. a small spreading of the Bloch angle k. In order 
to prepare phase space localized states, an optical lat- 
tice is carefully turned on. When the tunnelling exper- 
iment starts, the atomic density and the interaction be- 
tween atoms is sufficiently small, and the experiment can 
be analyzed as the interaction of individual independent 
atoms with the laser beams, i.e. using Hamiltonian (|^). 
However, the cloud of atoms remains cold enough, at a 
sub-recoil temperature, to prevent a large thermal broad- 
ening of momentum distribution that would destroy the 
signal. Because they start from very low temperature, 
these preparation techniques based on condensate ma- 
nipulation seem to allow a wider room to manoeuvre 
than those working with thermal clouds only. Adiabatic 
switching of the light potentials is not required and one 
can actually work with values of the classical parame- 
ters 7 and 9 which are far from the perturbative regime 
of an integrable system. 

By diagonalizing the evolution operator corresponding 
to eq. (H) over one period, we are not only able to repro- 
duce the oscillatory behavior of (p(t)) (see Fig. |^ (a)), 
but also we can study the spoiling effect of the thermal 
dispersion Ap oc -^/temperature and predict the maxi- 
mum allowed temperature (see Fig. || (b) and (c)). 1 If 
a denotes the width of the momentum distribution in re- 
coil momentum units, it can be shown §6. a] that it 
corresponds to a statistical mixture of Bloch states with 
Afc = a/2. Fig. U (a) corresponds to the ideal situa- 
tion where all atoms are prepared with a <C 1 about the 
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FIG. 3: Numerical simulation of the quantum evolution, in 
the conditions of the NlST experiment, i.e. 6 = 1.72 ,7 = 0.96 
and ftcft = 0.8 (compare to Fig. 4 (a) of jl3|). Starting 
at t = with a Gaussian wave packet whose Husimi func- 
tion is localized in one stable island (with a vanishing average 
momentum), we follow the average momentum (p) as time 
evolves. The stroboscopic measurements at times ir + 2mn 
(resp. 3n + 2mn) with m £ {0, 1, . . . ,40}, are plotted with 
the white (resp. black) circles. The tunnelling oscillations are 
clearly visible; the tunnelling period can be extracted from 
the typical time scale of the envelop: it is about 10 modula- 
tion periods. In the upper plot (a), we assume that a single 
Bloch angle k = is initially prepared (which implies a perfect 
phase coherence of the wavefunction across the optical lat- 
tice). The (thermal) dispersion of the Bloch angle washes out 
the signal: in case (b) , we take a momentum distribution with 
width Ap — a = 2Ak — 0.2 and in case (c) Ap = 0.4. In the 
latter case, the amplitude of the envelopes are so weak that 
this corresponds to an upper bound in temperature (about 
l/5th of the recoil temperature) at which tunnelling can be 
measured. 



k = Q. When a small but non vanishing a is introduced, 
some states of the quasi-energy bands with non vanish- 
ing k get involved and blur the tunnelling oscillations. 
For a — 0.2, the oscillation amplitude is reduced by a 
factor 2 and for a — 0.4 have nearly disappeared. There- 
fore, in this experiment, having a sub-recoil atom cloud 
is required. 

In the following we want to focus on tunnelling only 
and we will implicitly keep k = 0. 



B. Austin experiments [16| 



After the first version of the present p ape r was written, we learnt 
that an independent numerical work tej had obtained the same 
results as the ones in our Fig. ^ (a). But [0 do not consider 
the thermal effects and work always within theFloquet theory at 
k = instead of the Floquet-Bloch theory. In the present work 
we clearly demonstrate that thermal effects cannot be neglected 
and must be studied carefully when experiments are discussed. 



For a better understanding of the dynamics, one must 
go beyond the two-level model involving the symmetric 
and the antisymmetric states only. Other states must be 
taken into account and their influence can be enlightened 
when a classical parameter (7 or 9) or the quantum one 
(h Q g) is continuously varied. Two (quasi-)energies may 
exactly get degenerate if they belong to distinct symme- 
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try classes. If not, they may follow a so-called avoided 
crossing whose size reflects the direct coupling between 
the two states (more precisely the off-diagonal matrix el- 
ement of the coupling perturbation), but also the indirect 
coupling with other states. One of the keys of the chaotic 
tunnelling problems is to identify clearly the qualitative 
nature and the quantitative influence of indirect coupling. 
This is the background of Austin experiments. 

A third level is involved in a non negligible indirect 
coupling when its quasi-energy approaches the tunnelling 
doublet energies. This can be understood from perturba- 
tion theory as the leading term of the indirect coupling 
is proportionnal to the inverse of the energy difference. 
The two generic scenarii of the crossing of the doublet 
by a third state are shown in Fig. ^. Aside from the 




FIG. 4: When a third level is crossing a tunnelling doublet 
when a parameter of the Hamiltonian is varied, there is an 
avoided crossing between the third level and the member of 
the doublet with the same symmetry, while the other member 
of the doublet (with opposite symmetry) ignores the third 
level. Two generic scenarii exist: in case (a), the tunnelling 
splitting increases in the vicinity of the avoided crossing; in 
case (b), it decreases and vanishes at a specific value of the 
parameter. 

unavoidable ambiguous definition of the splitting, it ap- 
pears clearly that case a) corresponds to an increase of 
the tunnelling splitting during the crossing while, con- 
versely, case b) can lead to arbitrary small splittings since 
an exact degeneracy occurs. Thus, such a crossing by a 
third state produces a sharp variation of the tunnelling 
period that can be measured experimentally. This is ac- 
tually what is observed in the Austin experiments and 
can be confirmed by numerical experiments as shown in 
Fig. H When looking at the Husimi representation of 
the states, Fig. |[ one can immediately distinguish be- 
tween the tunnelling doublet and the third state suffi- 
ciently far from the crossing. As expected, the tunnelling 
doublet has Husimi representations localized in the sta- 
ble islands though they also spread in the chaotic sea. 
On a classical Poincare surface of section, it is easy to 
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FIG. 5: Numerical results obtained with the parameters of 
the Austin experiment, that is Hamiltonian (0), with 6 = 1 
and h a f{ = 0.33. The lower plot shows a part of the quasi- 
energy spectrum when 7 is varied. Thick lines show the two 
quasi-energies whose difference is the tunnelling splitting plot- 
ted in the upper plot. The two states are selected to have the 
largest localization of the Husimi function at the center of 
the stable islands. When a third level couples to the state 
that belongs to the same symmetry class, an avoided crossing 
can be seen and the definition of the doublet becomes neces- 
sarily ambiguous. There, some discontinuity in the selected 
state (and in the slope of the tunnelling frequency) cannot be 
avoided. The tunnelling splitting is shown as a function of 
7 in the upper plot, and is compared with the experimental 
results from ref. jl^] . The agreement is very good, which vali- 
dates the numerical approach. Around 7 = 0.2, a discrepancy 
is visible. This is precisely the "ambiguous" region where the 
dynamics cannot be reduced to a simple tunnelling oscillation, 
but at least three levels must be taken into account, leading 
to several relevant energy splittings. 



make the difference between chaotic and regular motion; 
how to transpose this distinction at the quantum level is 
not known with the large values of h c g used in both the 
Austin and the Nist experiments. Some classical struc- 
tures much smaller than the de Broglie wavelength are 
possibly present in some of the states 2 , but just looking 



2 From time to time, it is claimed |4o| that a simple matching 
between classical structures and quantum wave sub-Planckian 
structures was found but, as it was understood long time ago 
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NUMERICAL EVIDENCE OF A CHAOTIC 
TUNNELLING REGIME 




FIG. 6: (a), (b) and (c) show the density in gray scale of 
the Husimi functions associated with the three states which 
play an important role for tunnelling at 7 = 0.25 (parame- 
ters of the Austin experiment, as in fig. H. As expected, the 
two members of the doublet, (a) and (bj have their Husimi 
function localized about the two symmetric islands visible in 
the classical Poincare surface of section (d). (b) is strongly 
coupled to a third states shown in (c). h e g is too large to 
attribute any regular/chaotic character to the third state, (c) 
clearly plays a role in the enhancement of the tunnelling split- 
ting through an indirect coupling; this is thus a "assisted tun- 
nelling" mechanism, which cannot be unambiguously charac- 
terized as chaos assisted. 



It is one of the first success of quantum chaos to have 
shown that the energy levels of an integrable system ig- 
nore each others because they are localized on different 
classical tori, while in chaotic systems level repulsion is 
the rule. Following the discussion in the previous section, 
one may therefore expect that the average size of avoided 
crossings is increased when chaos is present. The fluctu- 
ations of a tunnelling process which should be narrow 
and sparse in an integrable regime should be broader, 
more numerous and possibly involving many states in a 
chaotic regime. We now illustrate this statement in the 
framework of the experimental atomic Hamiltonian (g) . 

There are two different ways of rendering chaos observ- 
able by quantum eyes: the first one consists in increasing 
the volume of the chaotic sea with the help of a classical 
parameter, the other one fixes the classical dynamics and 
decreases h c g. We will present both ways. 

CAT, RAT and all that... One may try to study the 
tunnelling fluctuations separately from the average be- 
havior. This average — in a somewhat vague sense — is 
increased because chaos diminishes the classical dynam- 
ical barrier 3 and this is the reason why the phenomenon 
can be called chaos assisted tunnelling. However, as far 
as only fluctuations are concerned (there can be an en- 
hancement or a decrease as well), the words "chaos as- 
sisted tunnelling" (CAT) [[To) may lead to confusion and 
we simply use "chaotic tunnelling". At last, Brodier, 
Schlagheck and Ullmo discovered what they called "res- 
onant assisted tunnelling" (RAT) 1 42 to describe an en- 
hancement of tunnelling due to an indirect process which 
involves one or several quasi-modes localized in the sec- 
ondary resonances surrounding the main symmetric is- 
lands. Up to now, RAT has been studied quantitatively 
in a purely quasi-integrable case but there are clues that 
it can be extended far beyond KAM theory. 



at the Husimi representation of the third state for h c ff ~ 1 
does not make it possible to attribute any chaotic or reg- 
ular character to it. It is only for much smaller h e g that 
chaotic or regular wavefunctions make sense. This is not 
surprising: the dichotomy between regular and irregular 
motion is a classical one and, at present, it can be ex- 
trapolated at the quantum level within the semiclassical 
regime only. Anyway, one must keep in mind that tun- 
nelling does only make sense in the semiclassical regime 
as well. One of the great merit of Austin experiments 
is to show for the first time a quantum tunnelling effect 
where an indirect process is involved, but we fill exager- 
ated to attribute any chaotic origin to it. 



by |4l| , the latter are generally washed out as soon as one try to 
measure global averages. 



A. 7 change 

Let us first take /i e ff ten times smaller than in the Nist 
and Austin experiments. It is much easier to do it nu- 
merically than experimentally, as it requires to increase 
the modulation frequency of the laser beams by one or- 
der of magnitude. We then follow the quantum states 
through the classical bifurcation shown in picture || and 
discussed in section IV A . Again, in order to calculate the 



tunnelling splitting Ae, we select the states that have the 



3 The question of defining an integrable system with respect to 
which the chaotic one should be compared is extremely difficult 
because of the exponential sensitivity of tunnelling on classical 
parameters. It is a much more serious problem than the way of 
defining an averaging procedure. 
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largest Husimi functions inside the islands. 
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FIG. 7: The tunnelling splitting as a function of 7, for 6 — 
1.724137 and h cB = 0.079638, that is 10 times smaller than 
in the nist experiments. Two regimes are clearly separated 
by the critical value 7 C ~ 0.625; these two quantum regimes 
correspond to two different classical regimes in the Poincare 
surfaces of section in Fig. |^. The smooth average decrease 
with sparse and narrow fluctuations (7 < y c ) corresponds 
to the case where the symmetric classical tori belong to the 
same regular island. The Hamiltonian can there be approx- 
imated by a simple integrable Hamiltonian, using the nor- 
mal form described in the appendix. The dotted line corre- 
sponds to the splitting calculated with this normal form and 
is in good agreement with the numerical result. In the sec- 
ond regime (7 > j c ), there are huge quantum fluctuations 
of the tunnelling splitting (and a slightly increased average 
value). Classically, this corresponds to tunnelling between 
unconnected symmetric islands, separated by a chaotic sea. 

Figure shows that, after a smooth decrease of Ae up 
to 7 = 7 C ~ 0.625 there is an abrupt change of regime. 
First, the mean value of Ae increases and second, many 
fluctuations appear which modify the splitting by several 
orders of magnitude. It is remarkable that this change of 
regime can be matched on the Poincare surfaces of sec- 
tion. The smooth tunnelling regime occurs the resonant 
tori belong to one quasi-integrable island. For 7 < ~y c , we 
are able to reproduce the main features of tunnelling by 
using an integrable approximation (see appendix). The 
decrease of tunnelling is directly understood in terms of 
the lengthening of the dynamical barrier. j c corresponds 
exactly to the point where the islands get disconnected. 
In addition, one can follow the bifurcation on the Husimi 
representations of the tunnelling doublet. 

One can hardly detect by eye any regularity in the 
chaotic regime but four large spikes in the range 7 C < 




FIG. 8: The tunnelling splitting as a function of h e g, for fixed 
values of the classical parameters, i.e. fixed classical dynamics 
(9 = 1.724137 and four values of 7). In the regular regime, 
an exponential decrease, as described by eq. ([!]) is visible, as 
a straight line with a negative slope in a logarithmic scale. 
The regime of fluctuations can be seen for 7 > 7 C when ft e ff is 
small enough for the de Broglie wavelength to be comparable 
to the size of the chaotic sea between the two islands. 



7 < 0.7. This can be traced back to the crossing by the 
same third state whose quasi-energy line is folded four 
times in the Floquet zone centered on the doublet. 



B. h c fi change 

Figure |§] shows Ae as a function of l/H e g for four values 
of 7. Here again, two types of regimes, a quasi-regular 
and a chaotic one, can clearly be distinguished. For the 
values of 7 where some substantial chaos is present in 
between the islands, by decreasing h c g, one gets to the 
chaotic regime. For 7 < 7 C , i.e. in the quasi-regular 
regime, one may enter in a chaotic regime but for far 
much low value of h e g since the chaotic layers between 
KAM tori are so thin that they are not even resolved in 
the Poincare surface of section given in the figures. 

On this plots, a purely exponential law given by eq. (Q) 
would produce a straight line with negative slope. It 
gives rather poor predictions in the chaotic regime and 
should be corrected even in the quasi-integrable regime 
in order to reproduce the fluctuations. Whenever such 
a fluctuation is due to a crossing by a quasi-mode (and 
not a state delocalized in the surrounding chaotic sea) , it 
might be reproduced by the resonance assisted tunnelling 
techniques. 
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VI. STATISTICAL SIGNATURE OF CHAOTIC 
TUNNELLING ? 

In order to reproduce quantitatively the statistics of 
the tunnelling splitting fluctuations in the chaotic regime, 
Leyvraz and Ullmo |l3| have introduced a random matrix 
model. The Hamiltonian can be split in two uncoupled 
components associated with the two even and odd sym- 
metry subspaces. The corresponding matrices are writ- 
ten as: 



jjeven 



(4 


v l "2 






vl 




vl 


Hi 


V 





(4) 



and 



H° 



dd 



Vl_ 



(5) 



where represent the energies of the doublet, Hi the 
Hamiltonian in the chaotic sea (modelled by a random 
Gaussian matrix) and v the indirect coupling. 

Neglecting direct tunnelling consists in taking Cq = £q ■ 
The central hypothesis is to consider all the v's as in- 
dependent variables with the same Gaussian distribu- 
tion. This is quite natural to treat all the other states 
on the same footing, as they are assumed to be chaotic 
states randomly delocalizcd in the chaotic sea. With 
these assumptions, the splitting distribution can be cal- 
culated and is given by a (truncated) Cauchy distribu- 
tion, see j|3|. For Hamiltonian (^), in each chaotic case 
where it has been tested, the Leyvraz-Ullmo prediction is 
in agreement with the numerical results (see jL4|). More 
surprisingly, we have found that the Leyvraz-Ullmo law 
gives correct predictions even when the classical dynam- 
ics is quasi-integrable (see Fig. |9| and the corresponding 
Poincare surface of section in Fig. [j]) and the de Broglie 
wavelength much larger than the chaotic layers. When 
looking at the splitting, one observes numerous and large 
fluctuations (over several orders of magnitude) which 
were supposed to characterize the chaotic regime. Con- 
versely, we have not been able to find a regular regime of 
fluctuations with a classically chaotic dynamics. There- 
fore one must conclude that even though two different 
regimes of tunnelling fluctuations can be identified un- 
ambiguously, classical chaos appears to be a sufficient 
but not a necessary condition for having numerous and 
large fluctuations governed by the Leyvraz-Ullmo law. 

These unexpected results are not explained at the 
present state of the theoretical approaches of chaotic tun- 
nelling. If it appeared, within future numerical or real 
experiments, that these results are not due to the pecu- 
liar properties of our system, it would definitely mean 
that theoretical studies should be all the more needed. 
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FIG. 9: Statistical distribution of the tunnelling splittings 
as ft e ff is varied in the quasi-integrable case corresponding to 
Fig. (b) {0 = 1, 7 = 0.018). Surprisingly, the numerically 
observed distribution (small circles) is in good agreement with 
the Leyvraz-Ullmo law (solid line) which is supposed to be 
valid in the chaotic it treats all the states coupled 

to the tunnelling doublet on the same footing. This clearly 
indicates that a Leyvraz-Ullmo distribution is not sufficient 
to characterize chaotic tunnelling. 



VII. CONCLUSION 

In this paper we have studied in great details the 
transition between a regular and a chaotic regime of 
tunnelling within a classical configuration that can be 
achieved experimentally. We have shown why h g has to 
be small enough if one wants to reach for the first time 
the chaotic regime in real systems. In recent experiments 
with cold atoms, it requires to increase the modulation 
period up to one order of magnitude (that is at least 
the MHz). 

Theoretically, there is also a lot of work to do if we 
want to understand and therefore predict the fluctua- 
tions quantitatively. The semiclassical regime requires to 
study carefully the dynamics with complex coordinates. 
However, the present work has shown clearly that the 
abrupt transition between a regular and a chaotic regime 
of tunnelling corresponds to a classical transition that can 
be identified very precisely. Tunnelling being a relevant 
concept in a semiclassical regime only, it is therefore not 
surprising that the future investigations on chaotic tun- 
nelling will have to keep track of the classical dynamics 
in one way or another. 
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APPENDIX A: INTEGRAL APPROXIMATION 
NEAR THE PITCHFORK BIFURCATION 

When one classical parameter is smoothly varied in 
a non-integrable Hamiltonian system, its periodic or- 
bits follow infinite fractal-like cascade of bifurcations 
through which they cannot be followed smoothly. How- 
ever, all the bifurcations can be classified according 
to a simple set of scenarios: for a given bifurcation 
of a given periodic orbit, the phase-space dynamics of 
the original parameter-dependent Hamiltonian can be 
uniformly approximated by an integrable parameter- 
dependent Hamiltonian which retains the relevant fea- 
tures only. Of course, because one cannot get rid of 
the chaotic dynamics, this approximation does make 
sense only locally, that is near the periodic orbit, in the 
neighborhood of the bifurcation. It is the object of the 
Hamiltonian normal form theory to classify the bifurca- 
tions and obtain the simplest form of the approximated 
integrable Hamiltonians (the so-called normal forms). 
Generically, i.e. when no constraint or symmetry is 
present, the one parameter Hamiltonian normal forms 
have been completely classified by Meyer |}0|, [ll], 
However, the bifurcation shown in Fig. || is out of the 
scope of Meyer's classification precisely because the time 
reversal symmetry plays a key role. Even if the Poincare 
surface of section near the origin (see Fig. || (b) cannot 
be distinguished from the one corresponding to Meyer's 
transitional bifurcation (see for instance Fig. 9 (b) 
in f3(j|l ), it is crucial to note that, in our case, two dis- 
tinct 27r-periodic orbits have emerged from the origin. 
In a transitional bifurcation, the two stable islands would 
correspond to the same 47r-periodic orbit. Therefore they 
would be classically connected to each other and would 
be irrelevant for tunnelling. Let us sketch briefly how to 
obtain the Hamiltonian normal form in our case. 

1. The first step is to find the value 70 of 7 at which 
the bifurcation occurs (9 is kept fixed). Such a bifurca- 
tion occur when the trace of the monodromy matrix M 7 
at the origin after one period is 2 and corresponds in 
the (— 27, 47#)-plane to the border of the even Arnold 
tongues (see Fig. 20.1 in |l6]]) (the transitional case oc- 
curs when trA/ 7 = — 2 at the border of the odd Arnold 
tongues). More precisely, the case shown in Fig. [2] cor- 
responds to 9 = 1.724137 (the experimental value in 
the Nist experiment) and the bifurcation takes place at 
70 — 0.564673. In the following, we use e = 7 — 70. 

Let us denote by y E (z) [resp. Z £ (z)] the solution of the 
Mathieu equation 

y"{x) + (47I? + 4 7 cos(2x))y(x) = (Al) 



such that y £ (Q) = and 3^(0) = 1 [resp. Z e (0) = 1 and 
Z' £ (Q) — 0]. The prime stands for the derivative with 
respect to x. For s = 0, it is straightforward to show 
that 



1 



(A2) 



where Z^(tt) ~ 1.480919 for 9 = 1.724137. 

2. The second step is to make a (linear) 2-7r-periodic 
canonical change of coordinate that eliminates the time- 
dependence in the quadratic part of Hamiltonian (|^), 
near the origin and uniformly in e. We are then led to 
the Hamiltonian 



-Sep 2 

2 p (A3) 



+higher order 27r-periodic terms 

where 01,(3,5 are e-independent coefficients. Only /3 = 
^d e y e (n) evaluated at e — will be relevant since a and 
6 can be eliminated by a suitable canonical change of 
coordinates following a method explained in J30| , section 
4.2]. In our case (3 ~ 2.008/-7T. We are therefore led to 
the following normal form of the quadratic part of the 
Hamiltonian: 



) 1 2 + 



1 1 

24^" 



+higher order 27r-periodic terms 



(A4) 



3. Following the same reasoning that lead to the tran- 
sitional normal form ||(], sections and 4.2], all higher 
order 27r-periodic terms, but the resonant terms of the 
form hk(s)p k , can be cancelled by a suitable canonical 
change of coordinates. Because of the time reversal sym- 
metry, the coefficient /13 vanishes identically and there- 
fore the leading order normal form is 



1 1 

247 



\(3ep 2 - \h 4 (0)p A - (A5) 



The explicit calculation of ^.4(0) is tedious but it can 
be estimated numerically by fitting the coordinate q = 
0,p = ± \Wk~fhi of the satellite 27r-periodic orbits for e > 
0. We obtain h A ~ 0.0320. 

It is far from obvious that the quantization of the nor- 
mal form (A5) will give a good approximation of the 
quasi-energies of the tunnelling doublet. Some discrepan- 
cies may rise from the fact that quantum physics is invari- 
ant under canonical transformations only at the leading 
order in h. When we swap p wit h q a nd change the sign of 
the energies, the normal form (A5) leads to a standard 
ID double-well problem whose quantum spectrum can 
be found by numerically diagonalizing the Hamiltonian 
written in a harmonic basis. The tunnelling splitting of 
the ground doublet is given by the dotted line in Fig. |?] 
and it agrees reasonnably well with the exact numerical 
result for 7 near 70. 
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